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ABSTRACT 



Context. It has been recently discovered that Cepheids harbor circumstellar envelopes (CSEs). RS Pup is the Cepheid that presents 
the most prominent circumstellar envelope known, the origin of which is not yet understood. 

Aims. Our purpose is to estimate the flux contribution of the CSE around RS Pup at the one arcsecond scale (~ 2000 AU) and to 
investigate its geometry, especially regarding asymmetries, to constrain its physical properties. 

Methods. We obtained near-infrared images in two narrow band filters centered on A = 1.644 and 2. 180/im (NB_1.64 and IB_2. 18, 
respectively) that comprise two recombination lines of hydrogen: the 12-4 and 7-4 (Brackett y) transitions, respectively. We used 
NACO's cube mode observations in order to improve the angular resolution with the shift-and-add technique, and to qualitatively 
study the symmetry of the spatially extended emission from the CSE with a statistical study of the speckle noise. 
Results. We probably detect at a 2cr level an extended emission with a relative flux (compared with the star in the same filter) of 
38 ± 17% in the NB.1.64 filter and 24 ± 11% in the IB.2.18 filter. This emission is centered on RS Pup and does not present any 
detectable asymmetry. We attribute the detected emission to the likely presence of an hydrogen envelope surrounding the star. 
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1. Introduction 

For almost a century, Cepheid stars have been used as distance 
indicators thanks to the well-known period-luminosity relation 
(P-L), also called Leavitt's law. A good calibration of this rela- 
tion is necessary to obtain an unbiased estimate of the distance. 
After the detection of circumstellar envelopes ( CSEs) around 
many Cepheids in the near- and mid-infrared (IKervella et alj 
I2006t iMerand et alJl2006l l2007t iKervella et alJl2009alT it was 
shown that these CSEs could bias the angu lar diameter determi - 
nation using near-infrared interferometry dMerand et al.ll2006l) . 
A circumstellar envelope could have an impact on the distance 
estimate if we do not take into account its flux contribution. It 
could lead to a bias both in interferometric and photometric mea- 
surements and therefore influence the accuracy of the estimated 
distance. Because Cepheids are used as standard candles to esti- 
mate distances in the universe, it is important to understand the 
properties of the CSEs around this kind of stars. 

We report new observations with VLT/NACO in the IB_2.18 
and NB_1.64 filters using a fast imaging mode in order to char- 
acterize the extended emission around RS Puppis (HD 68860) at 
sub-arcsecond angular resolution (»0.01 pc at RS Pup distance). 
Data and reduction methods are described in Sect. [2] In Sects. 
[3] and |4] we present two different data processing analyses: the 
shift-and-add method, and a method based on a statistical analy- 



sis of the speckle noise to enable the detection of the CSE in the 
speckle cloud surrounding the star's image. 

This paper is focuse d on RS Pup's envelope, which ha s been 
known for many years dWesterlundl 1 1 96 it iHavlenl 1 1 9721) for it 
is very extended; it has a spatial extension of around 2' in V 
band, which corresponds t o about one parsec at RS Pup's dis- 
tance (about 2 kiloparsecs: IKervella et al]l2008l) . Its origin is not 
precisely known; it is especially unclear whether it is made of 
material ejected by the star or of a material that predates the for- 
mation of the star. Its actual shape has been wi dely discussed: 



discre te nebular knots have been invoked by IKervella et al 



(2008) to estimate the distance using light echoes, whereas lFeast 



(120081) argued that a disk-like CSE will bias this kind of es- 
timate. Alterna t ely, a bipolar nebula has been considered by 
iBond & Sparks] d2009l) . making a ratio between two images of 
different epochs. We were able to spatially resolve the CSE close 
to the star. The morphology of the CSE is discussed in Sect. 14. 31 
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2. Observations and data reduction 

RS Pup was observed with the NACO instrument installed at the 
Nasmyth B focus of UT4 of ESO VLT. NA CO is an adaptive op- 
tics system (NAOS, iRousset et al. 200 3h and a high-resolution 
near IR camera (CONICA. iLenzen et al.ll200l . working as im- 
ager or as spectrograph in the range 1-5 /im. We used the S13 
camera for our RS Pup observations (FoV of 13.5"xl3.5" and 
13.26 + 0.03 mas/pixel. iMasciadri et al.ll2003l) with the narrow 
band filter NB_1 .64 (Aq - 1 .644 /mi with a width AA = 0.01 8 yum) 
and the intermediate band filter IB_2.18 (Aq = 2.180//m with a 
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width AA = 0.060 /mi). We chose the cube mode in order to apply 
the shift-and-add technique. We chose windows with 512x514 
pixels and 256x258 pixels respectively for the narrow and the 
intermediate band filters. Each cube (10 for RS Pup and 4 for 
the reference star) contains 460 images at 1.64 /mi and 2 000 
at 2.18/im. Each frame has the minimum integration time and 
readout type for this mode, i.e 109 ms and 39 ms respectively for 
the 1.64 and 2. 18 /mi cubes. With these DITs and filters we are 
in a q uasi-monochro matic regime with a very good speckle con- 
trast <lRoddieilll98lh . The large number of exposures allows us 
to select the frames with the brightest image quality. 

Data were obtained on 2009 January 7, with the atmospheric 
conditions presented in Table Q] The PSF (point spread func- 
tion) calibrator stars were observed just before or immediately 
after the scientific targets with the same instrumental configura- 
tions. Each raw image was processed in a standard way using 
the Yoriclfl language: bias subtraction, flat-field and bad pixel 
corrections. The negligible sky background was not subtracted. 

To compare the AO system behavior and artefacts on 
RS Pup, we used the star Achernar and its PSF reference star 
observed the same night. These data were also acquired in cube 
mode with the same camera in the NB_1.64 and NB_2.17 (Aq — 
2. 166 /mi with a width AA = 0.023 /mi) filters. Only one cube 
was acquired per filter, containing 20 000 frames for Achernar 
and 8 000 for the calibrator. The window detector size is 64x66 
pixels in both filters with an integration time of 7.2 ms. The 
slightly different second filter and the shorter exposure time 
compared with RS Pup is not critical for our analysis because 
we are interested in flux ratios. The smaller window size cuts 
the AO PSF wings, but we will show in Sect. |3]that this does not 
affect our conclusion significantly. 

Once the pre-processing was over, we carried out a precen- 
tering as well as a sorting according to the maximum intensity of 
the central peak to reject the images for which the correction was 
the least efficient. On the one hand, we applied the shift-and-add 
technique to obtain the best angular resolution possible. On the 
other hand, we applied a statistical study aimed at extracting a 
possible diffuse (i.e. completely resolved) component inside of 
the speckle cloud of the central core (the star). These methods 
are explained in Sects. [3]and|4] 

3. Shift-and-add process 

In this section we use the shift-and-add method and then com- 
pare the RS Pup images with another star observed on the same 
night for which we do not expect an envelope. 

3.1. RS Pup 

Both stars, RS Pup and HD 74417, are unresolved by the tele- 
scope: diffraction limits for an 8.2m UT are 55 mas and 41 mas 
in filters IB_2.18 and NB_1.64 respectively, while angular di- 
ameters are approximately one milli-arcsecond (1 mas) for both 
stars. 

The shift- a nd-add method was originally proposed by 
iBates & Cadvl (119801) as a technique for speckle imaging. 
Applied to NACO's cube mode, it enhances the Strehl ratio 
by selecting the frames that are the least altered by the atmo- 
spheric turbulence. It also reduces the halo contribution when 
used with adaptive optics by selecting the best AO-corrected 
frames. The efficiency of th e cube mode vs. the standard long ex- 
posure mode is discussed in lKervella et al.l (l2009bl) . [Girard et aD 



(1201 Oh showed for the RS Pup's reference star at 2. 18 /mi that 
by using this method we have an absolute gain of 9 % in Strehl 
ratio. 

Our processing steps are as follows: we select the 10 % best 
frames according to the brightest pixel (as a tracer of the Strehl 
ratio) in our 10 data cubes (4 for the calibrator). We then spa- 
tially resample them by a factor 4 using a cubic spline interpo- 
lation and co-align them using a Gaussian fitting on the central 
core, at a precision lev el of a few milliarcse conds (this method 
is described in detail in lKervella et al.ll2009ah . Each cube is then 
averaged to obtain ten final average images (4 for the calibra- 
tor). We finally compute the mean of the 10 average images of 
RS Pup (or the 4 for the calibrator). 

The result of this method is shown in Fig.Q]for the two filters. 
The two upper frames are HD 74417 and RS Pup as seen through 
the IB_2.18 filter, while the two lower are in the NB_1.64 filter. 
The spatial scale in all images is 3.3 mas/pix and we chose a log- 
arithmic intensity contrast with adjusted lower and higher cut-off 
to show details in the halo. On the 1.64 /mi image of RS Pup 
compared with the HD 74417 image we can see that there is a 
resolved circumstellar emission up to ~ 1" from the Cepheid. 
This feature is fainter in the 2. 18 /mi filter. In both filters, be- 
cause of this emission, diffraction rings are less visible around 
RS Pup than around HD 74417. 

For each final co-added frame (i.e the image shown in Fig.Q]) 
we computed a ring median and normalized to unity. The result- 
ing curves are presented in Fig.|2]for two filters. The error bars 
are the standard deviation of values inside rings and are plotted 
only for few radii for clarity. We notice a significant difference 
between the two curves, more pronounced in the NB_1.64 filter. 

3.2. Achernar 

We carried out the same processing for Achernar and its refer- 
ence star 6 Phe (HD 9362). These two stars are unresolved by the 
UTs (angular diameters are approximately 2 mas). With this win- 
dow size (64x66 pixels) the field of view is four times smaller 
than in the RS Pup's frames, and even if we can also see some 
Airy rings in all the images, the background is dominated by the 
AO halo. 

We plot in Fig. [2] the radial profiles of the stacked images, 
normalized to 10 (for clarity in the graph). 

3.3. Halo fitting 

We now investigate the presence of a circumstellar envelope by 
studying the radial profile of the stacked images we obtained 
for both RS Pup and its PSF calibrator and compare them with 
the observations of Achernar. The main idea behind this study 
is to compare the variations of encircled energy (the ratio of the 
core's flux to the total flux) to the expected variations caused by 
atmospheric condition changes. If RS Pup has indeed a circum- 
stellar envelope, its encircled energy will be significantly lower 
than that of our reference stars. 

We use an analytical function to fit the residual light out of 
the coheren t core (i.e . the h alo). We used a turbulence degraded 
profile from lRoddieJ(ll98lh 



4aloOO = / 



0.488 



6 \p 



(1) 



http://www.maumae.net/yorick/doc/index.php 



where the fitting parameters are the FWHM p and / proportional 
to the flux of the envelope. In this case, the parameter p corre- 
sponds to the equivalent seeing disk for long exposure. NACO's 
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Date 




UTC (start-end) 




Filter 


Name 


Seein g D IMM 


seeing,! 


AM 


Nbr 


r (cm) 


t (ms) 


2009-01 


■07 


00:30:41-00:31:56 


- 


NB.2.17 


HD 9362 


0.72" 


0.55" 


1.15 


1 


9.2 


3.1 


2009-01 


-07 


00:35:33-00:36:47 


- 


NB.1.64 


HD 9362 


0.77" 


0.62" 


1.15 


1 


9.0 


2.9 


2009-01 


-07 


00:51:16-00:54:16 




NB.2.17 


Achernar 


0.84" 


0.64" 


1.25 


1 


9.1 


2.3 


2009-01 


-07 


01:01:30-01:04:28 




NB.1.64 


Achernar 


0.83" 


0.63" 


1.27 


1 


9.1 


2.3 


2009-01 


-07 


04:47:17-05:00:46 


0.035 


IB^.18 


RSPup 


0.71 ±0.05" 


0.54 ± 0.04" 


1.03 


10 


13.6 ±0.4 


3.6 ±0.2 


2009-01 


-07 


05:01:16-05:13:17 


0.035 


NB.1.64 


RSPup 


0.68 ± 0.04" 


0.54 ±0.03" 


1.03 


10 


13.0 ± 0.1 


3.4 ±0.1 


2009-01 


-07 


05:21:03-05:26:27 




IB^.18 


HD 74417 


0.55 ±0.01" 


0.42 ±0.01" 


1.05 


4 


19.3 ±0.8 


5.1 ±0.4 


2009-01 


-07 


05:26:57-05:31:43 




NB.1.64 


HD 74417 


0.54 ±0.01" 


0.43 ±0.01" 


1.04 


4 


20.3 ±0.1 


5.1 ±0.1 



Table 1. Log of the NACO observations, seeingpjjjyjjyj are the seeing measurements in V from the Differential Image Motion Monitor 
(DIMM) station, while seeing^ are interpolated to our wavelength. The uncertainties denote the standard deviations over the number 
of cubes. The start-end UTC ind icated denotes t he tim e when the first and the last cube was recorded. The phase <p of the Cepheid 
was computed using (P, To) from lKervella et all (120081) . AM denotes the airmass, while Nbr is the number of cubes recorded, ro and 
to are the Fried parameter and the coherence time given by the AO real time calculator. 




HD 74417 1.64um 





Fig. 1. Final frames from shift-and-add process. In the two upper panels the PSF calibrator star HD 74417 and RS Pup are shown in 
the IB_2. 18 filter. The bottom panels show the same stars in the NB_1.64 filter. The scale intensity is logarithmic. 



resulting core, even for an unresolved source, is fairly complex 
and suffers from residual aberrations, this is why we decided to 
fit the halo at a large radial distance of the core. We will later 
define the core as the residual of the fit of the halo function to 
our data. 

We show the resulting fits in Fig. [4] They were applied for 
r > 0.26"at 2.18pm and 2.17pm while the range was for r > 
0.2"at 1.64pm for RS Pup (and its calibrator) and r > 0.26"for 



Achernar (and its calibrator). The black lines are our data (ring- 
median), while the gray lines represent the adjusted curves. The 
vertical scale is logarithmic and normalized to the total flux. We 
also evaluated the normalized encircled energy (£ e ncX i.e. the ra- 
tio of the flux in the core (i.e within 1.22/l/D) to the total flux. 
We can already see from these curves that the flux in the halo 
is different in the images of RS Pup and Achernar. These differ- 
ences can be linked to the emission of the CSE around RS Pup 
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Angular distance (arcsec) 




Angular distance (arcsec) 



Fig. 2. On the left are presented the radial profiles at 2. 18 pm for RS Pup, Achernar and their reference stars. On the right the same at 
1.64/im. Radial profiles were computed using the median value on one-pixel thick circular rings. The gray color denotes Achernar 
and 6 Phe profiles, shifted by a factor 10 for clarity. The scale is logarithmic for both axes. 




Fig. 3. Histogram of the encircled energy computed for 10% 
best frames. The solid lines represent a Gaussian fit. 



or a modification of the seeing between the two observations. 
This point will be discussed bellow with a statistical study of the 
encircled energy. 

For RS Pup in the NB_1.64 filter we found from our fit 
comparing the two FWHMs, i.e (p RS p up - Phd74417)/Phd74417, a 
change of halo size of 1 1 % and encircled energies £' enc (RSPup) 
= 0.327 and £ enc (HD74417) = 0.452 (see Table 0. This yields 
an encircled energy ratio s ~ 12%. The local DIMM station 
registered a seeing variation of around 26%. Different values 
can be explained by the station location, which is not inside 
the telescope dome, but a few tens of meters below, rendering 
it more suscepti ble to ground layer effects that are known to af- 
fect the seeing (iSarazin et al.ll2008l) . As a result the UT seeing 
is consistently better than the one measured by the DIMM. For 
IB_2.18 we show the fitting curves in Fig. [4] and the resulting fit 



parameters in Table [2] Our fit gives a seeing variation of around 
9 %, while the DIMM measurements give a 29 % variation, lead- 
ing to the same conclusion as before about the accuracy of the 
DIMM station seeing estimates. The encircled energies mea- 
sured in IB_2.18 are £ enc (RSPup) = 0.417 and £ enc (HD74417) 
= 0.517, giving an encircled energy ratio s ~ 80 %. 

We evaluated the uncertainties according to the variations of 
the encircled energy in the cubes. Because we already selected 
the 10 % best frames, we kept these samples and plotted the his- 
togram of the encircled energy to check the distribution function 
(Fig. 0. At the first order, the E em distribution can be approxi- 
mated by a Gaussian distribution. We then computed the relative 
standard deviation o-(E enc )/E enc for both stars. We assessed the 
uncertainty of the encircled energy ratio s with 



die) 



0~(£enc) 



RSPup 



C-(£enc) 



(2) 



HD74417 



We assume that the relative standard deviation of the 
encircled energy is a good approximation of the see- 
ing variations. We measured a relative standard deviation 
[cr(£ enc )/!e^] RS p U p ~ 8% at 2.18pm and -10% at 1.64/im 
while [cr(£ , enc )/£ , enc ] HD7 44 17 ~4% and ~6% respectively. We 
also evaluated the statistical uncertainties with a bootstrap 
method that gave cr ~ 0.5 % for both stars. These statistical 
uncertainties do not take into account the seeing variations. We 
therefore chose Eq. [2] as a reliable estimate of the final uncer- 
tainty. This gives cr(s)/s ~ 9% and ~ 12% respectively in 
IB _2. 18 and NB.1.64. 

For Achernar, the fit gives a seeing variation ~ 3% both 
in NB_1.64 and IB_2.17. The encircled energies measured in 
IB_2.18 are £ enc (Achernar) = 0.495 and £ enc (rJPhe) = 0.505, 
giving an encircled energy ratio s ~ 98 % (see Table |2|. In 
NB_1.64 we have 0.374 and 0.376 respectively, which gives 
s ~ 99%. The encircled energies are nearly identical, giving 
ratios close to 100%. Using Eq. [2] as previously, we found a 
relative standard deviation ~ 15 % in IB _2. 17 and ~ 11 % in 
NB_1.64. Because Achernar does not have an extended circum- 
stellar component in NB_1.64 and IB_2.17, we attribute these 
variations to the Strehl variations and the AO correction quality. 
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0.0 0.5 1.0 0.0 0.1 0.2 0.3 0.4 

r (arcsec) r {arcsec) 



Fig. 4. Halo fit for all filters with logarithmic ordinates. The fluxes are normalized to the total flux. Left panel: halo fit from the radial 
profile of images in the two filters for RS Pup. Right panel: halo fit from the radial profile of images in the two filters for Achernar. 
The fit in IB_2.18 was for r > 0.26" and in NB_1.64 for r > 0.2" for RS Pup and for r > 0.26" for Achernar. The dashed lines 
denote the reference stars, while the thick lines denote the Cepheid and Achernar. The fitted curves are plotted in gray. 



The encircled energy ratio seems stable between the two stars at 
the 2% level in IB.2.18 (see Table©. 

We can notice that the normalized curves of Achernar and 
its calibrator almost overlap, and this is what is expected for un- 
resolved stars without envelope taken in same atmospheric and 
instrumental conditions. It is not the case for RS Pup (see Fig. [2] 
and |U). The additional emission coming from the envelope of 
the Cepheid adds a component that is clearly not negligible. In 
the same atmospheric conditions and AO correction between the 
star and its calibrator, the encircled energy ratio should be close 
to 100 %. That is what we obtain from Achernar's data. Scaling 
RS Pup's images to the same window size (i.e. r < 0.42") 
slightly changes £ enc because we normalized to the total flux, 
but we still have a discrepancy of 100% for the encircled en- 
ergy ratio (s ~ 79 % and s ~ 87 % respectively in NB_1 .64 
and IB_2.18). A smaller window on RS Pup makes the total flux 
lower because we do not integrate all the envelope flux, conse- 
quently the normalized encircled energy is larger and therefore 
the E enc ratio as well. 

The drop in E enc of RS Pup compared with the statistics over 
Achernar's data and the fact that this drop is larger than the one 
directly estimated from the seeing halo variations, motivate the 
following section where this drop is analyzed under the assump- 
tion that it is caused by CSE emission. 



3.4. Encircled energy variation 

In this part we evaluate the flux of the envelope using the fact that 
the theoretical encircled energy ratio should be 100% and that an 
additional component to the flux should decrease this ratio. 



filter 


name 


P(") 




/ 


F 

^enc 


NB.1.64 


RSPup 


0.46 ± 0.01 


0.58 


±0.01 


0.327 




HD 74417 


0.43 ± 0.01 


0.43 


+ 0.01 


0.452 




Achernar 


0.49 ± 0.01 


1.00 


+ 0.01 


0.374 




<5Phe 


0.46 ± 0.01 


0.97 


±0.01 


0.376 


IB_2.18 


RSPup 


0.56 ± 0.01 


0.45 


±0.01 


0.417 




HD 74417 


0.50 ± 0.01 


0.34 


±0.01 


0.517 


NB.2.17 


Achernar 


0.52 ± 0.01 


0.59 


±0.01 


0.495 




rSPhe 


0.54 ± 0.01 


0.61 


±0.01 


0.505 



Table 2. Best-fit results for the residual halo where, p is the 
FWHM and / is a parameter proportional to the flux of the halo 
normalized to the total flux. E eDC represents the flux in the central 
core (i.e. within 1.22A/D) to the total flux. Note: The fitted /halo 
analytical function is not normalized to /, i.e J Ih R i rdr + f. 



The encircled energy of RS Pup can be written as 



where F is the total flux of the star, without the extended flux, 
F CO h is the flux in the coherent core, F env is the emission of the 
environment and a is the fraction of the envelope flux that lays 
in the core. For simplicity, we will assume a « 1, which means 
that we assume that the CSE is significantly larger than the core. 
Accordingly we can write 

£ _ F coh /F 

enc * i + F env /r 
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Assuming that the encircled energy is stable in NACO for 
our observing conditions, we can calibrate the previous formula 
to obtain the extended emission flux: 

£enc(Cepheid) 1 

s = ~ 

E en c (reference) 1 + F(env)/F(Cepheid) 

The relative emission /^(envVfXCepheid) can be estimated 
with the encircled energy ratio from our measured values 
(Table |2|. We can also take the uncertainty values assessed in 
the previous section, defined as being the relative standard de- 
viation on 10% of best frames in the cubes. We found a relative 
extended emission around RS Pup contributing to 38 + 17 % of 
the star's flux in the NB_1.64band and 24 + 11 % in IB_2.18. 

4. Morphological analysis 

We now present a qualitative analysis of the envelope morphol- 
ogy. In this analysis, based on the statistical properties of the 
noise owing to speckle boiling, we presented the ratio cr(F)/F as 
a function of the factor (g). This parameter is intrinsic to speckle 
temporal and spatial structures, is not expected to depend on az- 
imuth, and is proportional to the envelope flux. We mapped this 
factor, and any azimuthal asymmetry should be attributed to the 
morphology of the envelope. 

4. 1 . Analysis method 

The principle of the method is to recover the flux of the enve- 
lope in the short exposures. The difference between a simple un- 
resolved star and a star surrounded by an envelope is that the 
envelope will prevent the flux between the speckles from reach- 
ing zero. 

Statistical descriptions of stell ar speckles have 
been proposed by man y authors dRacine et al 
iFitzgerald & Graham! |2006|) . who showed that speckle noise 
is not negligible. Th e distribution of sp e ckle s tatistics in AO 
were also studied b y ICanales & Cagigal (Il999h . For example, 
iRacine et al.l d 1999b (hereafter Ra99) showed that the speckle 
noise dominates other noise sources within the halo of bright 
stars corrected with AO systems. We know that the total 
variance in the cube short exposures after the basic corrections 
(i.e corrected for bias, bad pixels, flat-field and recentered) is 

cr 2 {r) = oi(r) + C7^ h (r) + cr 2 Rm + oj behi (r), (3) 

where cr 2 is the variance of the speckle noise, cr 2 h the variance 

of the photon noise, cr^ 0N the variance of the readout noise and 

°"strehi tne n °i se owing to Strehl ratio variations. We neglected 
the variance of the sky photon noise, which is negligible in our 
wavelength for bright targets. Because the Strehl ratio variations 
occur in the coherent core and because we are particularly inter- 
ested in the halo part, we omit these variations in the rest of this 
section. 

We present for one calibrating cube the variance observed 
in the top panels of Fig. [5] for both filters. We plotted the radial 
profile computed using median rings on all azimuths. The dashed 
curve denotes the variance of the readout noise. We can estimate 
some of the components in Eq. [3] from our data. For example, 
from the stacked image that we obtained by summing the cube, 
we obtain the core and halo profile, as we did in Sect. [3] We then 
obtain the photon noise variance using its Poissonian property, 
knowing that the variance is equal to the mean. We define for the 



stacked image the radial intensity profile I(r) and we estimate the 
term <x 2 h of Eq. [3] with 

where N is the number of frames in the cube (1999 frames in 
IB_2.18 and 457 in NB_1.64) and y the gain of the detector in 
order to convert noises in electron unit. The readout noise is es- 
timated from the outer part of the images (see in Fig. [5] the top 
panels). We estimate the speckle noise from the total variance of 
the cubes as 

a 2 (r) = cr 2 (r) - (r 2 ph (r) - cr| 0N , 

In the bottom panels of Fig. [5] we present the relative im- 
portance of the different noise sources. We see that the different 
variances show a different radial behavior, and we see in par- 
ticular that the photon noise of the core never dominates. More 
precisely, we can define two regimes: 

- cr(r) ^ <x s (r) for r < r\, in this regime, the speckle noise 
dominates the other variances (cr s > 3 ct RO n) 

- oif) - ^oi(r) + o-2 QN for r > r\ 

withrj = 0.22" in NB.1.64 band and n =* 0.24" in IB_2.18. The 
fairly weak (oc A~ l2 l 5 ) wavelength dependency of the speckle 
noise may explain the same order of magnitude between the limit 
radius in IB_2.18 and NB_1.64. 

Indeed there is another regime in the range < r < r\ where 
the Strehl variations occur. The PSF core is mainly affected by 
these variations. As said previously, we are interested in the PSF 
halo to detect the envelope between the speckles, which is why 
we neglect this term. 

4.2. Speckle noise invariant 

We use an interesting property of the speckle noise, namely 
that its variance is proportional to the squared flux in the cube 
(Ra99), F/cr is therefore an invariant. More precisely: 

cr s = g(r,A,At,S,...)F±(A), 

where the function g is a function depending on the wavelength 
A, the integration time At, the Strehl ratio S and atmospheric 
parameters. 

We checked this property by plotting in Fig. [6] the mean g 
function for each star, i.e we evaluated a g function for each 
cube that we then averaged (the error bars correspond to the stan- 
dard deviation between the cubes). We took care to remove the 
flux contribution from the circumstellar envelope with the values 
found in Sect. [3] Obviously there is no clear difference between 
the two stars except in the core. This can be explained by the 
Strehl ratio variations we omitted before in the core. We can as- 
sess these variations by evaluating the encircled energy for each 
frame of each calibrating cube (here we take 100% of frames, 
without selection) and by using it as an estimator of the Strehl ra- 
tio in function of time. By fitting a Gaussian distribution to E enc 
for each star and both filters, we es tima te the relative standard 
deviation for RS Pup to be cr{E enc )l E enc ~ 15 % in IB_2.18 and 
~ 17 % in NB.1.64. For HD 74417 we found ~ 8 % and ~ 9 % 
respectively at 2. 18 /mi and 1.64/iin. 

The total invariance of the g function, i.e for the complete 
radius, depends on the Strehl ratio variations. Nevertheless, be- 
cause these variations occur in the coherent core, we can neglect 
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Fig. 5. Left panel: at the top we plot the variance (circles) and the average value (triangle) of one data cube of HD 74417 at 2. 18//m. 
Bottom: contribution of the different noise sources. Right panel: the same for the 1.64/im filter. The scale for all plots is logarithmic 
in ordinates. The dotted curves show the 3cr limits of the readout noise, while the gray curves denote the halo fit. 
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Fig. 6. Mean g function for each star. The error bars correspond to the RMS. The bottom plot represents the absolute relative 
difference. We can see a larger difference in the core owing to the Strehl variations. From r > r\ the functions are similar at a 9 % 
level in NB.1.64 and at 6 % in IB_2.18. 



them if we go far away from the central part. For r > r\ we 
found a mean absolute relative difference (see Fig. [6]l of 5 % in 
IB.2.18 and 8 % in NB.1.64. So for r > r { the invariance of the 
g function is verified with a good accuracy. 

We can measure g using our reference star. Assuming it is 
stable in time, we can define the observable T as the average 
of the flux of a given pixel in the cube, divided by the standard 
deviation. Defining F cll \ = cvFrs Pup = aF* and making the dif- 
ference between r sc i and F ca i (where the index sci and cal refer 
to RS Pup and HD 74417 respectively) we have 



■r C2 



aF+ 



cal 



(4) 



+ cr; 



RON 



All variables in Eq. [4] are functions of the pixel coordinates 
(x, y), but we omit this notation for clarity. 

The parameter a can be estimated with the results from 
Sect. [3] The flux ratio between the two final images of RS Pup 
and HD 74417 is a function of a and F(env)/F*: 



F(RS Pup) 
F(HD 74417) 



F* + F a 
a F+ 



1 + F em /F* 
a 



(5) 



By measuring the total flux in our final images (Sect. [3]) and 
using F env /Fi, estimated in Sect. l3.4l we get a = 1.06 ± 0.13 at 
1.64 pm and a = 0.96 + 0.08 at 2.18 ptn. We also used the total 
flux variations in our 10 (for RS Pup) and 4 (for the reference) 
averaged images to assess the uncertainties, but they turn out to 
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be small (< 2 %) compared to F em /Fi, uncertainties. As a first 
approximation we can simplify Eq.[4]with a ~ 1 to obtain 

g= ^ :■ (6) 

^8 2 *Fl + 4qn 

£ is thus directly proportional to the flux of the CSE. We 
expect this function to be null if there is no envelope. The de- 
nominator is a radial parameter, therefore any spatial asymmetry 
in £ will be linked to an asymmetry of the CSE. 

From our different data cubes we computed a mean £ by sub- 
tracting to each r sc j a mean r ca i, i.e: 

%(x,y) = - V(r sci j - T ca i) and F ca i = — V r calii , 

where n and m are the number of cube for RS Pup and HD 74417 
respectively. This process provides us with an average image 
proportional to the flux of the CSE in each filter. 

4.3. Symmetry of the CSE 

In Fig. |7]we show the f function for the NB.1.64 and IB_2.18 
band. For radii larger than r\ , we expect the function to be null if 
there is no CSE (see Eq.[6]i, whereas we clearly obtain a positive 
image, indicative of a likely detection of a CSE. The central part 
of the image with r < r\ (white circle for r - r\) denotes the part 
where the previously defined function g is not invariant because 
of Strehl ratio fluctuations. It is difficult to estimate and therefore 
does not supply reliable informations close to the star. 

This method is useful for studying the morphology of the 
CSE, in particular whether or not it has a central symmetry, like 
one expects for a uniform shell structure or a face-on disk; or 
an asymmetrical shape, like in an inclined disk-like structure for 
instance. Despite adaptive optics artefacts, the images in Fig. [7] 
seem to present a uniform intensity distribution in agreement 
with a shell or a face-on disk structure. More precisely, we esti- 
mated a symmetry level using a residual map obtained by sub- 
tracting 90 degree rotated versions of themselves from the im- 
ages presented in Fig. 7. The result of this subtraction is shown 
in Fig. [8] In the initial image we assessed the average value in 



a ring, while in the subtracted image (0° - 90°) we compute an 
average value in a small window (see Fig. [8]). Making the ratio 
of these two averages (i.e. the average value in the small window 
over the average value in the ring), we exclude a departure from 
symmetry larger than 7 % in NB_1.64 and 6 % in IB.2.18. 




Fig. 8. Estimate of the symmetry level. Top panels: on the left 
is plotted the £ parameter at 2.18//m with a ring of radius 
r\ < r < ri — 0.92" where we computed a mean value. On the 
right we plot the difference of the initial (0 °) and 90 ° rotation 
images where we cut a small window to assess a mean residual 
value. The ratio between the mean residual value and the mean 
value give us a level of symmetry. Bottom panels: same process 
at 1.64yum with a ring of radius r\ < r < r% = \.\" . The scale 
intensity is logarithmic. 
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5. Conclusion 

We present a study of the close environment (1000-2000 AU) of 
the Cepheid RS Pup using near-infrared AO-assisted imaging. 
We used two different techniques to 1) measure the photometric 
emission due to the CSE and 2) obtain an image proportional 
to the CSE through a statistical analysis of the large number of 
short exposures in which speckles are still visible. 

Using a simple image modeling method, based on the un- 
resolved PSF calibrator, we estimated (Sect. [3} a likely pho- 
tometric emission of 38 ± 17 % and 24 ± 1 1 % of the star flux 
in the NB_1.64 and IB_2.18 bands respectively for RS Pup's 
CSE. Apart from the actual interest to better understand the mass 
loss in Cepheids, the presence of such an extended emission 
may have an impact on the distance estimate with the Baade- 
Wesselink method. 

Using an original statistical study of our short exposure data 
cubes, we qualitatively showed (Sect. |4| that the envelope in 
the two bands (NB_1.64 and IB_2.18) has a centro-symmetrical 
shape (excluding a departure from symmetry larger than 7 %), 
suggesting that the CSE is either a face-on disk or a uniform 
shell. 

Our observations were obtained in two filters isolating the 
Brackett 12-4 (1.641yum, NB.1.64 filter) and Brackett y 1-A 
(2.166yt/m, IB_2.18 filter) recombination lines of hydrogen. We 
therefore propose that our observations show the presence of hy- 
drogen in the CSE of RS Pup. The star is naturally relatively faint 
in these spectral regions, because of deep hydrogen absorption 
lines in its spectrum. This appears as a natural explanation to the 
high relative contribution from the CSE compared to the stellar 
flux. Another contribution to the CSE flux may also come from 
free-free continuum emission, or Rayleigh scattering by circum- 
stellar dust. However, this contribution is probably limited to a 
few percents of the stellar flux. The high relative brightness of 
the detected CSE therefore implies that the spectral energy dis- 
tribution of the detected CSE is essentially made of emission 
in hydrogen spectral lines instead of in the form of a contin- 
uum. RS Pup's infrared excess in t he broad NB_1 .64 and IB_2. 1 8 
bands is limited to a few percents dKervella et al.ll20~09ah . which 
excludes a continuum CSE emission of 40% in the NB_1 .64 filter 
and 24% in the IB_2.18 filter as measured in our NACO images. 

Our discovery is also remarkably consistent with the non- 
pulsating hydrogen absorp tion componen t detec ted in the visi- 
ble Ha line of RS Pup bv iNardetto et al.l (120081) . together with 
a phase-dependent component in emission. The extended emis- 
sion we observe in the infrared is therefore probably caused by 
the same hydrogen envelope as the one that creates the Ha emis- 
sion observed in the visible. It should als o be mentioned that our 
result is consistent with the proposal by iKervella et al.l (l2009al) 
that the envelope of RS Pup is made of two distinct components: 
a compact gaseous envelope, which is probably at the origin of 
the observed NACO emissi on, and a very large (» 2') and cold 
(w 45 K) dusty envelope dKervella et al]l2008l) that is mostly 
transparent at near-infrared wavelengths owing to the reduced 
scattering efficiency. These two envelopes appear to be largely 
unrelated because of th eir very differen t spatia l scales. However, 
it should be noted that iMarengo et al.l d2009l) recently detected 
large envelopes around several Cepheids using Spitzer. The pres- 
ence of hydrogen close to the Cepheids of their sample could be 
investigated at high angular resolution using the observing tech- 
nique presented here. 
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